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Abstract 

We study the one-dimensional Riemann problem for a hyperbolic system of three con¬ 
servation laws of Temple class. The system is a simplification of a recently propose system 
of five conservations laws by Bouchut and Boyaval that models viscoelastic fluids. An im¬ 
portant issue is that the considered 3x3 system is such that every characteristic field is 
linearly degenerate. We study the Riemann problem for this system and under suitable gen¬ 
eralized Rankine-Hugoniot relation and entropy condition, both existence and uniqueness of 
particular delta-shock type solutions are established. 


1 Introduction 

In this work, we show existence and uniqueness of delta shock wave for the Suliciu relaxation 
system mm 

Pt + ( pu)x = 0, 

< (pu)t + ( pu 2 + s 2 v) x = 0, (1) 

k (pv)t + (puv + u) x = 0, 

with s = const. > 0. The concept of delta shock wave is a generalization of an ordinry shock 
wave. This generalization was introdduced by Korchinski in the year of 1977 in his PhD thesis 
H3j. He considered the Riemann problem for system 

f«t + (i« 2 ) x = 0 , 

\vt + {h uv ) x = °- 
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Motivated by some numerical results, he constructed the unique Riemann solution using gen¬ 
eralized delta functions to obtain singular shocks satisfying (J2J) in the sense of distributions. 
After that, in 1994, Tan, Zhang and Zheng established in [20 ] the existence, uniqueness and 
stability of delta shock waves for a viscous perturbation of system ([ 2 ]) where the first equation 
is replaced by ut + ( u 2 ) = £u xx , with e > 0. Other works dealing with delta-shock waves are 
due to Ercole m who in 2000 obtained a delta shock solution as a limit of smooth solutions 
by the vanishing viscosity method. Sheng and Zhang m discussed the Riemann problem for 
pressureless gas system. In 2005, Brenier [2] considered the Riemann problem for the Chaplygin 
gas system. Other works for the Chaplygin gas system can be found in mmm ■ Some more 
recent works on delta shocks for general hyperbolic conservation laws are due to Danilov and 
Mitrovic Ei, where they described delta shock wave generation from continuous initial data by 
using smooth approximations in the weak sense. The works mentioned above on the delta shock 
solutions are for particular cases of 2 x 2 systems. For particular 3x3 non-strictly hyperbolic sys¬ 
tems with repeated eigenvalues, we mention the works of Panov and Shelkovich m, Shelkovich 
m and Cheng [5] . For a strictly hyperbolic system, see [9] for partial results of delta shock 
solution. For new developments of delta shock waves, its applications and historical notes is Ell- 

Coming back to system ©, we can deduce that the associated eigenvalues are given by, 
Ai = u— A2 = u and A3 = u+ ^ where the corresponding Riemann invariants are R\ = s 2 v—su, 
R 2 = v + i and R 3 = s 2 v + su. From the expressions for the eigenvalues and the Riemann in¬ 
variants we obtain Ai = ^ — SR 2 , X 2 = 55 (R .3 — R\) and A3 = SR 2 — We can see that 
system m is linearly degenerate. 


Recently, Lu et al. [15] showed existence of solutions for the Cauchy problem associated to 
the Suliciu relaxation system m with bounded initial data 


(p(0,x),u(0,x),u(0,x)) = (p 0 (x),u 0 (x),v 0 (x)), po(x)>0, 16 R, (3) 


subject to the following conditions: 

(HI) The functions po, uq and vo satisfy 

Cl < uq(x) — SVo(x) < c 2 , C 3 < u 0 (x) + sv 0 (x) < c 4 , 


v 0 (x) + 


1 

£ + p 0 (x) 


> c 5 , 


where e > 0 is a small and the constants c, , i = 1,..., 5, satisfy C 5 — C4 2s ci > 0. 

(H2) The total variations of ito(x) — svo(x) and uq(x) + svq(x) are bounded. 

The existence result for the Cauchy problem ([I])-([3]) includes solutions in vacuum regions jT5j, 
Theorem 1], When po(x ) > p > 0, the condition HI becomes: 

(HI) The functions po, uq and vo satisfy ci < uq(x) — sv o(x) < C2, C3 < uq(x) + suo(x) < c 4 , 
and uo(x) + > C5, where c,; , i = 1 ,..., 5 , are constants satisfying C5 — C4 2s Cl > 0 - 

In [9], we showed uniqueness of solutions for the Cauchy problem associated to the Suliciu 
relaxation system (E)-© with po(x) > p > 0. Moreover, Theorem 2 in [9] ensures that when 
we consider vq(x) = — , problem (E) is reduced to following Chaplygin gas system 

\p t + ( pu) x = 0 , 

\{pu)t + ( pu 2 - £) x = 0 . 







The classical Riemann problem associated to the Suliciu relaxation system (H|)-([3|) has been 
extensively studied, for instance in BEIHEI. In |10| . it is shown uniqueness for the generalized 
Riemann problem for the Suliciu relaxation system. 

In the present paper, we show the existence and uniqueness of solutions for the Riemann 
problem associated with the Suliciu relaxation system with initial data 


(p 0 (x),u 0 (x),v 0 (x)) 


( Pl,ui,vi ), if x < 0 , 
(p r , U r ,V r ), if X > 0, 


(4) 


in which the left and right constant states ( pi,ui,vi ) and (p r ,u r ,v r ), with pi,p r > 0 , satisfy the 
conditions H2 and Xi{pi,ui,vi) > A 3 (p r ,u r ,v r ), i.e., they do not satisfy condition HI globally. 
Under this assumptions we have the following situations: 

1. ( pi,ui,vi ) and (p r ,u r ,v r ) satisfy locally the condition HI, i.e., 

• ( Pi,ui,v t ) satisfy 


aq < ui — svi < «2) «3 < ui — svi < a4 and vi~\ -> a5 

Pi 

where a*, i = 1,..., 5, are suitable constants satisfying 0:5 — ° 4 ~ Q1 > 0, 

• ( p r ,u r ,v r ) satisfy 

/3i < u r — sv r < @21 P 3 < u r — sv r < 04 and v r H-> ^5 

Pr 

where /%, i = 1,..., 5, are suitable constants satisfying /3s — ■ 4 »yi > 0. 

• Let ci = min{ai,/3i} and C 4 = max{aq,/?i}. Then there are Ci, i = 1,... ,4 such that 


ci < 


r ui - svi 
— sv r 


< C2 


and 


(ui + svi 1 
\u r + sv r J 


< c 4 , 


but is not possible to find a constant C 5 such that C 5 — C4 2s Cl > 0 . 

2. Only (p r ,u r ,i y) satisfy locally the condition HI, i.e., ( pi,ui,vi ) satisfy 


/?i < ui ~ svi < f3 2 , Pz<ui- svi < 04 and vi H- > /3 5 

Pi 

where ft, i = 1,..., 5, are suitable constants satisfying ft — < 0. 

3. Only ( pi,ui,vi ) satisfy locally the condition HI, i.e., (p r ,u r ,v r ) satisfy 

ft ft u r — sv r < ft, ft < u r — sv r < ft and v r H-> ft 

Pr 

where ft, i = 1,..., 5, are suitable constants satisfying ft — < 0. 

4. Neither ( pi,ui,vi ) nor (p r ,u r ,v r ) satisfy the local condition HI. 







De la cruz et al. [9j studied only the first case. The idea used by the authors in [9] can be 
extended to the other cases. 

We considered all the cases for which existence and uniqueness of solutions can be guarantied 
for the Riemann problem associated to the Suliciu system. In |9j it is shown that if the initial data 
m, with po(x ) > p > 0, satisfies HI and H2, then, the Riemann problem has a unique solution. 
It is also shown that HI and H2 imply the Lax shock condition Xi(pi,ui,vi) < A 3 (pi,m,vi). 
Therefore, the Riemann problem for the Suliciu relaxation system has a classical solution in 
the region Ti = {(p,u,v) : p > 0 and Xi(pi,ui,vi) < A 3 (p,u,v)}. In this paper we show that 
the Suliciu relaxation system presents nonclassical solutions in the sense of delta shocks in the 
region 

T 2 = {(p,u,v) : p > 0 , Ai (pi,ui,vi) > A 3 (p,u,v) and (ui - u ) 2 > s 2 (vi - v)(l/p - 1 /p{j}. 

It remains to study if we can obtain some kind of solution in the region (0, 00 ) x R 2 — (TiU^). 
This object of our current research. 


2 Delta shock solutions 


Denote by BM(M) the space of bounded Borel measures on R. the definition of a measure 
solution of Suliciu relaxation system in BM( R) can be given as follows. 

Definition 2.1. A triple (p, u, v ) constitutes a measure solution to the Suliciu relaxation system, 
if it holds that 


a) p G L°°((0,oo),HM(R)) n C((0, 00 ), 1L~ S (R)), 

b) u € L°°(( 0 ,00), L°°(R)) n C(( 0 ,00), LT S (R)), 


c) ueL-((0,oo),L-(R))nC((0,oo),^ s (R)), s > 0, 


d) u and v are measurable with respect to p at almost for all t € ( 0 , 00 ), 
and 

h = / 0 °° / R Ot + u<t>x) dpdt = 0 , 

( h = / 0 °° f R u(<pt + u(f> x ) dpdt + f 0 °° f R s 2 vcj) x dxdt = 0, (5) 

„ 13=r Ir v ^ t + d P dt +/o°° Ir udxdt= 
for all test function <j) € C'^°(R + X R). 

Recall that a two-dimensional weighted delta function w(s)5l supported on a smooth curve L 
parameterized as t = t(s), x = x(s) (c < s < d) is defined by ( w{s)5l , (f>(t, x)) = w(s) 4 >(t(s),x(s)) ds 

for all 0 € Cq°(R 2 ). We next present the definition of a delta shock wave. 


Definition 2.2. A triple distribution ( p,u,v ) is called a delta shock wave if it is represented in 
the form 


(, p,u,v)(t,x ) 


(. Pi,ui,v t ){t,x ), x<x(t), 

< (w(t)S(x - x(t)),us(t),g(t)), x = x(t), 
(p r ,U r ,V r )(t,x), X>x(t), 


( 6 ) 


and satisfies Definition ED where ( pi,ui,vi)(t,x ) and (p r ,u r ,v r )(t,x) are piecewise smooth 
bounded solutions of the Suliciu relaxation system ©■ 




We set = us(t) since the concentration in p need to travel at the speed of discontinuity. 
Hence, we say that a delta shock wave © is a measure solution to the Suliciu relaxation system 
0 if and only if the following relation holds, 


£ ^ = Mt), 

T = ~[p]us(t) + [pu ], 

d w(t)u s {t) = - [pu}u s (t ) + [pu 2 + s 2 r>], 

k dw(t Jl j(t} = -[pv\us(t ) + [puv + u}. 


(7) 


In fact, for any test function cj) € C^°(M + X M), from ([5]), we obtain 

h = J j {fa + u(j) x ) dpdt = J |-ti, 5 (t)[p] + [pu] - dW ^ | (j)dt, 

h = J l^-us(t)[pu] + [pu 2 + s 2 u] - | ( ]f. an d 

h = J | -us (t) [pv] + [puv + u] - dU | 4>dt. 

Relations (0 are called the generalized Rankine-Hugoniot condition. It reflects the exact rela¬ 
tionship among the limit states on two sides of the discontinuity, the weight, propagation speed 
and the location of the discontinuity. In addition, to guarantee uniqueness, the delta shock wave 
should satisfy the admissibility (entropy) condition A 3 {p r ,u r ,v r ) < us(t) < Ai (pi,ui,vi). 


Now, the generalized Rankine-Hugoniot condition 0 is applied to the Riemann problem 0- 
© with left and right constant states U- = (p-,U-,V-) and U+ = (p + ,u + ,v + ), respectively, 
satisfying the condition H2, the fact A 3 (p+, u + , v + ) < Ai(/?_, it_, u_) and 


(u- — u + ) 2 > — (v- — v + ) —— (v- — v + ). 


( 8 ) 


P+ P- 

Therefore, the Riemann problem 0)-© is reduced to solving 0 with initial data 

t = 0, 2 ( 0 ) = 0, w{ 0) = 0, 5 ( 0 ) = 0, (9) 

under entropy condition 


s s 

u .(_ -|-< U$(t) < U- -. 

p+ p- 

From 0 and (I9|) j it follows that 

w(t) = ~[p\x(t) + [pu\t, 
w(t)us(t ) = — [pu]x(t) + [pu 2 + s 2 -y]t, and 
w(t)g(t) = ~[pv\x(t) + [puv + u]t. 


( 10 ) 


( 11 ) 


Multiplying the first equation in (flTT) by ug(t) and then subtracting it from the second one, we 
obtain that 

[p]x(t)us(t) — [pu]us(t)t — [pu\x{t) + [pu 2 + s 2 v]t = 0 , 


[p] 2 U \ r, [pu 2 + s 2 v\ 2 


Y® (*) - \pu]x(t)t + 


t z = 0 , 


that is, 


d 

dt 


2 











which is equivalent to 


[p\x 2 (t) — 2 [pu\x(t)t + [pv 2 + s 2 v]t 2 = 0. (12) 

From m, one can find us(t) := us is a constant and x(t ) = ust. Then, (fl2|) can be rewritten 

[p\u 2 - 2[pu]us + [pu 2 + s 2 v] = 0. (13) 

When [p\ = p_ — p + = 0, the situation is very simple and one can easily calculate the solution 


_ U-+U+ 2 M 

US ~ 2 ^ b 2 p_ [u] ’ 

x(t) = Ust, 

w(t ) = P-(U- — U + )t, 

m = 

which obviously satisfies the entropy condition cm From condition 

2 M / 1 


P- 


< -(Ar(Z7_) - A 3 {U + )Y < -M(Ar(C/_) - A 3 (£4)) 


and 


us - I U- 


P- 


u + + 


p- 


- \U- - 


p- 




p-m 


< o. 


Similarly, we can deduce that 


because 


P- 


us — \u+ H-= - \u- -- H-+ s 


P- 


P- 


p-m 


>o, 


2H / 1 


1 


s-j- < - X 3 (U + )Y < -M(Ar(t/_) - A 3 (U+)). 

When [p] = p_ — p + Y 0, the discriminant of the quadratic equation (fl3l) is 
A = 4 [pu] 2 — A[p][pu 2 + s 2 u] = p-p + [u] 2 — s 2 [/9 ][u] > 0 
and then we can find 

( _ WjPy/Ag 

Us ~ [p] 


(S*) = 


x(t) = 

w(t ) = ±y/A/At, 

, , _ - [pu] [pv}±[pv\ yj A/4 +[p] [puv+u] 

[p]\/A74 


(14) 


(15) 


Next, with the help of the entropy condition (1101) . we will choose the admissible solution 
from (USD. Observe that by the entropy condition and since the system is strictly hyperbolic, 
we have that 

s s s s 

u .|_- < u + < u .|_ H-< U- -< U- < U- H-. 

P+ P+ P- P- 













Observe that 


-[p\\i(p-,u-,v-) + [pu] = p+ 
-[p]^3(p+,«+,w+) + \pu\ = P- 



( s 
[ u+ -v + 

(u + + — 

V p+ 


> 0 , 

> 0 , 


[p](\i(p-,u-,v-)) 2 - 2[pu]Xi(p-,u-,v-) + [pu 2 ] + s 2 [u] = 


P+ \ u - — u+ 


2 s 2 


+ — + s 2 [u] < 0 , 
P- 


\p\(X 3 (p + ,u + ,v + )) 2 - 2[pu]X 3 {p + ,u + ,v +) + [pu 2 ] + s 2 [u] = 
P- (u- —u + -— ^ — — + s 2 [w] > 0, 

V P+J P+ 

then, for the solution given in (1151) by (S’ -1- ), we have 

us — Ai(/9 _, U-, V-) < 0 and u$ — X 3 (p + , u + , v + ) > 0, 


which imply that the entropy condition (flOl) is valid. When Ai(p_,u_,u_) = X 3 (p+,u + ,v + ), we 
have trivially that Xi(p-,U-,V-) = us = X 3 (p + ,u+,v+). 

Now, for the solution ( S~ ), when p_ < p + we have 


us - X 3 (p + ,u + ,v + ) 


P-(X 3 (U~) - X 3 (U+)) + V [puf - \p][pu 2 + s*v] 
\P] 


and when p_ > p + , that 


, P+ (Xi(U-) - Xi(u+)) + VM 2 -[d[p « 2 + * 2 «] . , 

us - Ai(p_,u_,u_) =-r-p-> ( 

\p\ 

showing that the solution ( S~) in (fl5l) does not satisfy the entropy condition (fTOl) . 
Thus we have proved the following result. 


Theorem 2.1. Given left and right constant states ( pi,ui,vi) and (p r ,u r ,v r ), respectively, such 
that satisfy the condition H2, X\(pi,ui,vi) > X 3 {p r ,u r ,v r ) and dSJ), that is, 

(Ui - u r ) 2 > S 2 (vi - v r )/p r - S 2 (u/ - V r )/pi- 


Then, the Riemann problem <m> o admits a unique entropy solution in the sense of measures. 
This solution is of the form 


( P,u,v)(t,x ) 


' (pi,ui,vi), 

< (w(t)S(x - u s t),u s ,g{t)), 

„ iPr i U r , V r ), 


if x < ust, 
if x = u 6 t , 
if x > ust, 


where us, w(t) and g(t ) are show in (fTTl) for [p] = 0 or ( S + ) in (|T5D for [p\ / 0 . 







3 Numerical illustration 


In this section, we show numerical evidence of delta shock solution for the Suliciu relaxation 
system using the Lax-Friedrichs method. In the numerical test, with s = 1, we consider the 
initial data given by 


(Po,uo,vq)(x) 


(9,5,14/5), if x < 0, 
(1,3,2), if x > 0, 


the spatial discretization parameter for N = 1780 points and a constant CFL = 0.1969889. 

The numerical results at final time t = 0.1 is displayed in Figure [l] The exact solution at time 
t is 


{p,u,v){t,x) 


(9,5,14/5), if x < ust, 

< (At5(x — ust),u$, Bt), if x = u$t, 
k (l,3,2), if x > ust, 


with us 


2175-737 A 
475 ’ A 


2737 

75 


and B 


75+29737 

10737 


4 Conclusions 

In this paper we obtain delta-shock waves for the Riemann problem associated to the Soliciu 
relaxation system. This work complements and extend the partial results concerning delta- 
shock waves for the Riemann problem for the Suliciu relaxation system that were obtained in 
M- We mention that in [HE] the classical Riemann problem for the Suliciu relaxation system 
was solved. In our analysis we consider all possible cases for the existence and uniqueness 
of delta-shock solution for the Suliciu relaxation system and by the Lax-Friedrichs method we 
show numerical evidence of such solutions. We identified particular regions where we can present 
classical solutions and delta shock waves. However, we observed that there are regions in which 
the Riemann problem for this system does not have classical solution neither delta shock waves, 
which encourages us to seek solutions in some different sense to be able to analyze this missing 
regions. 
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Figure 1: Numerical solution. Top: (a) The initial data po(x) = p(0, x). (b) Numerical solution 

of p at time t = 0.1. Middle: (a) The initial data uq(x) = u(0,x). (b) Numerical solution of u 

at time t = 0.1. Bottom: (a) The initial data vq{x) = u(0,x). (b) Numerical solution of v at 

time t = 0.1 



























